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Проведено моделирование системы неразрушающего контроля 
дефектов в твердых телах. Рассматриваются обратные 
геометрические задачи теории упругости для плоской 
прямоугольной области по реконструкции круговых полостей 
и трещин, выходящих на поверхность тела. Дополнительной 
информацией для решения этих задач является набор первых 
четырех собственных резонансных частот. Решение обратных 
задач основано на минимизации функционала невязки между 
измеренной входной информацией и рассчитанной в ходе 
численного решения прямых задач с заданными параметрами 
дефектов. В качестве инструмента решения прямых задач 
используется метод конечных элементов, реализованный в 
программе НехРПЕ. функционала 
осуществляется с помощью генетического алгоритма (ГА), 
реализованного в разработанной программе САЕЕММОТ. В 


Минимизация 


работе описан алгоритм этой программы и настройки ГА, 


которые используются в численных экспериментах. 
Приведены результаты этих экспериментов по определению 
параметров дефектов (координат центра, радиуса, координат 
выхода трещины на поверхность и ее размер). Эти результаты 
показывают достаточность дополнительной информации для 
преодоления некорректности задачи, а также высокую 
эффективность предложенного алгоритма, как в точности 


определения параметров дефектов, так и в скорости их поиска. 


Ключевые слова: генетический алгоритм, метод конечных 
элементов, неразрушающий контроль. 
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Введение. Неразрушающий метод контроля (НРК) — контроль свойств и характеристик объекта, при котором 
не должна быть нарушена пригодность объекта к использованию и эксплуатации. Способ НРК прочности заключается 
в том, что исследуемая конструкция или материал не подвергается механическим разрушениям, контроль 


осуществляется косвенно путем измерения и математического анализа физико-механических величин, отвечающих за 


Работа выполнена в рамках инициативной НИР. 
"Е-пай: зо]о\Яеуагс @ этаЙ.сот, тезоцие @ этаЙ.сот 
"Тре гезеагсь 15 допе т Фе #гате оЁ Фе шдерепаеп( В &). 


Вестник Донского государственного технического университета 2016, №2(85), 5-12 





прочностные свойства конструкции или материала. 

Существующие методы неразрушающего контроля, такие как радиоволновой, оптический, капиллярный, 
тепловой, а также радиационный хорошо зарекомендовали себя и успешно используются в производстве и не только. 
Но с развитием техники и планомерным увеличением удельной производительности ЭВМ возникают новые 
возможности для проведения неразрушающего контроля. Одним из таких способов является использование 
эволюционных алгоритмов при распознавании измеряемой акустической информации об инспектируемом объекте. 
Среди этих методов наибольшее распространение в задачах дефектоскопии получили применение искусственных 
нейронных сетей (ИНС) и генетических алгоритмов (ГА) на этапе минимизации функционалов невязки между 
измеренной и рассчитанной информацией. Так применению ИНС при идентификации трещиноподобных дефектов в 
различных конструкциях посвящены работы [1-4] и др. Разработку и использование ГА в задачах идентификации 
механических свойств упругих тел и реконструкции дефектов можно найти в работах [5-8], в которых в качестве 
измеренной информации используются данные о вибрации тел или их резонансных частотах. 

В настоящей работе показана работоспособность сочетания метода конечных элементов и генетического 
алгоритма для решения задач НРК, на примере определения параметров внутренней круговой полости и 
прямолинейной трещины выходящей на поверхность упругого прямоугольника. Идентификация дефектов основана на 
минимизации функционалов невязки между данными о наборе первых четырех его резонансных частот. 

Постановка обратной задачи. Рассмотрим установившиеся колебания с круговой частотой ‹ упругого тела, 


занимающего конечную односвязную область У с границей 5. Граница области может быть описана объединением 


непересекающихся областей 5 = 5, (155 (53, причем на 5, заданы кинематические, а на 5, — силовые граничные 
условия. Множество 5, — это неизвестные границы дефекта (трещина, выходящая на поверхность, внутренняя 
полость), свободные от напряжений. 
Амплитудные характеристики этих колебаний удовлетворяют системе дифференциальных уравнений 
линейной теории упругости [9]: 
= ре’и. О; = Суми 1 1=1,2,3, ХЕ\, (1) 
и граничным условиям прямой краевой задачи 


—. 20 — — — 
и; 5 =и; о 5 =бут, |, = р, Е 15, =0, (2) 


где и;, Е, — компоненты векторов смещений и напряжений; с„, — компоненты тензора упругих постоянных; р — 
ПЛОТНОСТЬ. 
В обратной задаче внутренние границы области 5. являются неизвестными, информацией для их определения 
является набор собственных частот колебаний: 
(2 = {61 ,0,..., 6, }. (3) 
На практике наиболее точно определяются набор первых частот, вместе с этим информативность этих данных 
существенно зависит от вида и размера дефекта, что может быть выражено зависимостью 0; = 0,(5.), [=1,2,...п. 
Определение внутренней границы 5. сводится к нахождению абсолютного (нулевого) максимума функционала, 
задающего умноженную на —1 норму (в некотором пространстве Г, ) разности между измеренными резонансными (3) и 
найденными из решения краевой задачи (1)—(2) собственными частотами ©; = ©; (5; ай Те о заданная 
конфигурация дефекта, при этом У — область изменения искомых параметров, характеризующих 
дефект 5., определяется их физическим смыслом. В настоящей работе в качестве такого функционала используется 


выражение «взвешенной» нормы разности, при которой влияние на его значение изменения всех используемых частот 


В определенной степени одинаково: 
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Е($;)=- ы ‚_ 9$) ‚ ЗЕМ, (4) 


21| 6; (53) 


при этом область \ изменения параметров, характеризующих область 5., известна и определяется размерами 


исследуемого образца. 

Описание программного обеспечения. Для решения обратных геометрических задач идентификации 
дефектов разработан и программно реализован комплекс Сепейс А!еогт \уи Ешце Нетепе Мефоа ш М№п- 
Резбисйуе Тезипе (САРЕММОТ). Данный программный продукт является связующим звеном между модулем ГА 
«АГ::Оепейс::Рто» разработанного Зи тееск1 ГлКаз7 [10] и конечно-элементном пакетом, причем в настоящей работе 
был использован НехРПЕ. Результатом работы САЕЕММПОТ является набор параметров, обеспечивающий экстремум 
функции многих переменных. 

Система написана на скриптовом языке программирования РЕВТ, [11]. Выбор в пользу скриптового языка был 
продиктован необходимостью оперативно менять целевую функцию, что является не тривиальной задачей для 
классических языков, таких как Оер1, С++ и т. д. САЕЕММОТ представляет собой модульный РЕВГ-скрипт. Целевая 
платформа для запуска приложения: Глпих. 

Входными данными для САЕЕММОТ является файл, содержащий в себе описание целевого функционала и 
предельные значения варьируемых переменных. 

Выходными данными модуля САЕЕММОТ является текстовый файл, содержащий в себе таблицу с 
результатами расчета. Решением является строка с максимальным значением целевой функции. Простота 
представления выходных данных позволяет производить визуальный и статистический анализ практически в любом 
программном продукте, ориентированном на обработку табличных данных, начиная от офисных пакетов Мисгозой 
ОЁйсе и Гл6геОЁйсе и заканчивая узкопрофильными пакетами обработки данных, таких как Рага\Мле\ и У1$П. 

Примеры решения задач неразрушающего контроля. 


Ниже приведены результаты численных экспериментов по идентификации круговой полости и 
прямолинейной трещины, выходящей на поверхность прямоугольника под прямым углом. Идентификации подлежат 
дискретные параметры, характеризующие координаты центра и радиус окружности, или абсцисса точки выхода 
трещины на поверхность и ее длина. В качестве значений, характеризующих точность идентификации искомых 
параметров, рассматривались относительные ошибки е ‚ которые вычислялись по формуле: 


и с: (5) 
А 
где А — целевое значение искомого параметра; В — варьируемое расчетное значение этого параметра. 
Первым шагом минимизации функционала (4) с помощью ГА является настройка его параметров. 
Проведенные численные эксперименты показали, что наиболее рациональными в смысле «точность — время» 


являются следующие параметры генетического алгоритма: 


Таблица 1 
Параметры генетического алгоритма 


128 
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Идентификация полости в прямоугольном теле. Описанную методику проиллюстрируем на задачах 
идентификации внутренних полостей в прямоугольнике, совершающим собственные колебания в рамках задачи 


плоской деформации. Размеры прямоугольника: длина 0,02 м, ширина 0,01 м, материал — сталь, модуль 
10 2 
Юнга ЕЁ =20х10` Н/м”, коэффициент Пуассонау = 0,29. Правая сторона прямоугольника жестко закреплена, 


граница внутренней полости свободна от напряжений; Хи У — координаты центра; К — радиус окружности. 
Решение обратной задачи идентификации кругового отверстия осуществляется на основе минимизации 
функционала (4). В качестве инструмента решения прямой задачи использовался метод конечных элементов (МКЭ), 
реализованный в пакете НехРОЕ. Роль «измеренных» резонансных частот выполняет набор первых четырех 
собственных частот, рассчитанных при заданных параметрах дефекта. 
В табл. 2 приведены результаты расчета с помощью САЕЕММОТ для 10 примеров прямоугольников с 
отверстиями, которые были выбраны случайным образом. 


Таблица 2 
Поиск внутренней полости 


ох [| о» | омю | ое = 

у [| 008 | оз | ол = 
ох [| 0095 | обо | 008 — 
у [| 0» | ом |0 = 
к [| ож | ом | 009 — 
ох [| 00% | 0075 | 00% — 
к [| 00з [| обжз | от = 
ох [| 000 [| обы | 00% — 
у [| юз | обюз | ом 
Г [| 0042 | ом» | 00% — 
ох [| 0055 | обяю | оз — 
к [| от [| ок | 049 — 

ох [| ом | обюм | оп = 

у [| 00455 [| 005 | 006 — 

к [| 905% [| 005 | ом — 

ох [| 005 [| обжим | 04% 

у [| 00 [| оз | 098 — 

к [| мя | омуя | ож — 

ох [| 0008 [| оз | 008 — 

у [| 008 [| 08 | оби — 

о [| оз [| обв | оби 

ох [| ож | обв | 00% — 

у [| 00655 [| обв | оби — 

к [| 09я | о | 00% — 

ох [| 002 | обе | оби = 

т у [| особи | обвя | 00% — 
к [| 900 [| 900% | 000 — 


м 


аа а во 5 | 89| < | >< [55| < | 54 | 25 
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В 0.0015 | Х 0.004 у 0.003 
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6) 0.0014 В) 0.0037 0 м 
0.0036 -0,001 
0.00135 
0.0035 -0.002 
|| 10 19 28 37 46 55 | 10 19 28 37 46 55 
0.0013 о 0.0034 > -0.003 


Рис.1. а) Целевая функция; 6) А — радиус окружности; в) иг) Х и У — координаты центра окружности 


На рис. 1, а представлены графики изменения целевого функционала РЁ (4) и параметров К, Х и У. На рис. 
1, 0, в иг проиллюстрированы зависимости от М числа эпох ГА для задачи | из таблицы 2. Анализ этих результатов 
показывает, что приемлемая точность идентификации достигается уже при 40 эпохах. 


Прямоугольное тело с трещиной, выходящей на поверхность. Вторым примером применения описанной 
выше методики является идентификация выходящей под прямым углом на поверхность прямолинейной трещины в 
прямоугольнике длиной 0,05 м, шириной 0,01 м, из того же, что и в п. 3.1, материала. Половина правой стороны 
прямоугольника закреплена, ХО — координата, где трещина выходит на поверхность, У/ — глубина трещины (рис. 2). 





ло (ХО; УГ 





Рис. 2. Модель тела с трещиной 


В качестве входных данных были использованы первые три собственные частоты колебаний контролируемого 
объекта. В табл. 3 представлены результаты идентификации 10 примеров трещин, выходящих на нижнюю сторону 
прямоугольника. 

Таблица 3 
Поиск трещины, выходящей на поверхность 


№ 


[61 № [| 00% [ 0095 | 0915 — 
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Рис. 3. Целевая поверхность 
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Рис. 4. а) Целевая функция; 6) ХО — координата, где трещина выходит на поверхность; 
в) У] — глубина трещины 

На рис. 3 представлен пример целевой поверхности функционала (4) (значения ХО, У1, Е нормированы и 
изменяются в пределах [0;1]). На рис. 4 проиллюстрировано поведение целевой функции и параметров в зависимости 
от числа эпох М№ для 1-ой задачи из таблицы 3. Быстрая сходимость алгоритма (20 эпох) обясняется структурой 
целевой поверхности. 

Выводы. Результаты исследования состоят в следующем: 

® Определены рациональные параметры настройки ГА, представленные в табл. 1. 

е Для определения трех параметров полости (п. 3.1) с достаточной степенью точности требуется не менее 
четырех собственных частот колебаний. При этом удовлетворительная точность идентификации достигается при 40 
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поколениях ГА. 

е В задаче идентификации трещины (п. 3.2) достаточной информацией оказался набор первых трех 
резонансных частот, при этом целевая поверхность имеет один экстремум, для нахождения которого достаточно 20 
поколений ГА. 
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